Structure, Stability, Edge States and Aromaticity of Graphene Ribbons 
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We determine the stability, the geometry, the electronic and magnetic structure of hydrogen- 
terminated graphene-nanoribbons edges as a function of the hydrogen content of the environment by 
means of density functional theory. Antiferromagnetic zigzag ribbons are stable only at extremely- 
low ultra-vacuum pressures. Under more standard conditions, the most stable structures are the 
mono- and di-hydrogenated armchair edges and a zigzag edge reconstruction with one di- and 
two mono-hydrogenated sites. At high hydrogen-concentration "bulk" graphene is not stable and 
spontaneously breaks to form ribbons, in analogy to the spontaneous breaking of graphene into 
small-width nanoribbons observed experimentally in solution. The stability and the existence of 
exotic edge electronic-states and/or magnetism is rationalized in terms of simple concepts from 
organic chemistry (Clar's rule). 

PACS numbers: 71.15.Mb, 71.20.Tx, 73.20.At 
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While two-dimensional graphene [lj exhibits fascinat- 
ing properties such as a relativistic massless-dispersion 
for its charge carriers and ballistic transport on large 
distances, its gapless spectrum makes it unsuitable for 
direct application as a channel in field-effect transis- 
tors and other semiconducting devices. This problem 
might be overcome, and thus open the way to a break- 
through carbon-based electronics, by designing graphene 
nanor ibbons (GNRs) in which the lateral quantum 
confinement opens an electronic gap, function of the 
ribbon width. Recent works 0, |3| have shown that 
the use of lithographic patterning of graphene samples 
can yield GNRs. However, the reported ribbons have 
large widths, between 15 and 100 nm, small electronic 
gap, up to 200 meV, and are characterized by a signif- 
icant egde-roughness. For practical applications, larger 
gaps, and correspondingly narrower and smoother-edged 
GNRs (width w < 10 nm), are desirable. An alterna- 
tive experimental route [H, @] , that makes use of chem- 
ical methods such as solution-dispersion and sonication, 
has shown that graphene sheets spontaneously break into 
ribbons of narrow width and smooth edges. Ideally, a 
combination of the two methods could lead to the design 
of GNRs having edges of controlled orientation, spatial 
localization, and electronic properties. 

In this regard, the knowledge of the structural and 
energetic properties of the possible edges, as well as of 
their thermodynamic stability, is crucial to achieve the 
experimental control necessary for technological applica- 
tions Q . A huge number of theoretical works has recently 
appeared in the literature on GNRs @, H @, 0, EH, E , 
13 . H ]. In those studies, ribbons are usually chosen with 
the two fundamental edge geometries, zigzag and arm- 
chair, and dangling bonds are saturated with a single 
hydrogen per carbon atom. Zigzag GNRs have caused a 
great resonance in the scientific community, due to their 
surprising electronic and magnetic properties. In fact, 
zigzag GNRs feature magnetism due to spin-polarized 



electronic states localized on the edge [7, 8] and they 
possibly turn to half-metal under high external electric 
fields [3, [l(|. All this has fired the hope to use GNRs 
in future spintronic-devices 0, 11 1 , even if the robustness 
of the spin polarization in presence of defects has been 
questioned [14] . Surprisingly, no attempt has been made 
so far to study the energetics and thermodynamics of rib- 
bon edges having different hydrogen terminations, as a 
function of the chemical potential typical of the experi- 
mental conditions. 

In this work we report density functional theory (DFT) 
calculations of the energetics and structure of various 
hydrogen-terminated graphene-nanoribbons edges as a 
function of the hydrogen content of the environment. We 
show that magnetic nanoribbons are stable only at ex- 
tremely low hydrogen-concentrations, challenging to ob- 
tain experimentally. They are very reactive and thus 
unlikely to be stable at ambient conditions. The most 
stable structures at reasonable hydrogen-concentrations 
are non-magnetic and have a very low chemical reactiv- 
ity. They are thus likely to be metastable well beyond 
their region of thermodynamical stability. We also show 
that at high hydrogen-concentrations graphene is not 
stable and spontaneously breaks to form hydrogenated 
edges. This is analogous to the spontaneous breaking of 
graphene into small-width nanoribbons observed experi- 
mentally in organic compounds solution [H, 0] . Finally, 
the stability and the existence of exotic electronic edge- 
states and/or magnetism are fully rationalized in terms 
of simple concepts from organic chemistry (Clar's rule). 

We studied GNRs with different edge-terminations us- 
ing periodic super-cell geometries with lattice-parameter 
L (periodicity) along the edge, the two edges having the 
same configuration. Details of the total-energy calcu- 
lations are in note [lj|. Zigzag edges are denoted as 
Znin 2 ...n x where n; stands for the number of hydrogen 
atoms on a given site, and X is the number of adja- 
cent edge-sites within the periodicity L. Armchair edges 
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FIG. 1: A: Scheme and electronic band-structure of the five most stable hydrogen-terminated edges of a graphene nanoribbon. 
Carbon-carbon bonds are represented with the standard notation, while hydrogen atoms are the filled circles. The structures 
are periodic along the ribbon-edge with periodicity L. The gray area corresponds to the electronic-bands allowed in "bulk" 
graphene. The dashed line is the Fermi level. B, C: Other stable armchair and zigzag terminations. D: Further examples of 
possible edges. E: Standard representation of the benzenoid aromatic carbon ring. 



are denoted as a mini2 (a mims , m3m4 ) for supercells contain- 
ing one (two) hexagon column(s), where nij indicates the 
number of hydrogen atoms bonded to the i th carbon site. 
Some examples are in Fig. [T] Other structures which do 
not fit this notation, z(57)y, z(600)yjy, and a(56)i, are 
also in Fig. [1] Again, the subscripts indicate the number 
of hydrogens on a given site. 

In Tab. U we report the computed zero-temperature 
edge formation-energy per length: 

^ 2 = ^(E ribb -NcE blk -^E H2 y (1) 

where E r , E blk and Eh 2 are the total energy of the 
ribbon super-cell, of one atom in "bulk" graphene and of 
the isolated H2 molecule [3]. Nc (Nh) are the number 
of carbon (hydrogen) atoms in the super-cell. £ h 2 can be 
used to determine the stability of different structures as a 
function of the experimental conditions (l7| . In presence 
of molecular H2 gas, at a given chemical potential pn 2 , 
the relative stability is obtained by comparing Gh 2 = 
£h 2 — phph 2 /2, where pn — Nh/(2L). At the absolute 
temperature T and for a partial H2 pressure P [HI, 

Ph 2 = H°{T) - H°(0) - TS°(T) + k B T\n (J^j , (2) 

where H° (S°) is the enthalpy (entropy) at the pres- 
sure P° = 1 bar obtained from Ref. [19]. In pres- 
ence of monoatomic-hydrogen gas one should use Gh = 
£h - PhPh, where £ H = £h 2 — Ph X 2.24 eV, from DFT. 

The most stable structures are thus obtained by com- 
paring Gh 2 or Gh- G is linear in p, the slope of the 
line being determined by ph- For a given value of p 
the stable structure is the one with the lowest value of 
G, thus, by increasing p (i.e. going to an environment 



TABLE I: Formation energy (Eq. [TJ and hydrogen density 
(pH = Nh/(2L)) for all the studied edges. 
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richer in hydrogen) the favorable structures will be those 
with higher hy drog en-density pu- This concept is visu- 
alized, as usual [17], by plotting G vs. p in Fig. [2j where 
we report the five most stable structures. In Fig. [31 we 
report the same plot for the most stable edges in the two 
distinct families of zigzag and armchair GNRs. 

Some important conclusions can be driven from Fig. [2] 
First, the magnetic zi edge, widely studied as a potential 
component for future spintronic devices [^, [ll[ , and the 
z (57)oo, proposed in [lj|, are stable only at extremely low 
hydrogen-concentrations. Second, at reasonable hydro- 
gen pressure (at ambient conditions, the partial H2 pres- 
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FIG. 2: Formation energies versus chemical potential for the 
five most stable edges. Vertical lines distinguish the stability 
regions. The alternative bottom (top) axes show the pressure, 
in bar, of molecular Hb (atomic H) gas corresponding to the 
chemical potentials at T = 100, 300, and 600 K. 
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FIG. 3: Formation energies versus chemical potential for the 
most stable armchair (top) and zigzag (bottom) edges. 



sure in air is ~ 5 x 10 -7 bar) the most stable edges are an , 
Z211, and a22- Finally, Fig. [2] shows that graphene is not 
stable at high hydrogen-concentration (G becomes neg- 
ative) and spontaneously breaks to form hydrogenated 
edges. Similarly, it has been experimentally observed 



that small-width nanoribbons spontaneously form in or- 
ganic compounds solutions @, Q . 

We indicate in Tab. Q] along with the energetics, 
whether the proposed edges are magnetic or not, i.e. 
whether a finite magnetization density is observed within 
DFT calculations. In Fig. [2l we also report the elec- 
tronic bands of the ribbons and we show the bands of 
two-dimensional "bulk" graphene projected on the one- 
dimensional Brillouin zone of the edge (gray area in 
Fig. [2]). As currently done in surface physics, we call 
"edge states" (ESs) those electronic states localized at 
edges, having an energy which is not allowed in the bulk 
(i.e. those bands which are outside the gray area in 
Fig. [2J. Such states are localized near the edge. More- 
over, in Tab. U we indicate as "metallic edges", those 
having partially- filled ESs, i.e. having ESs which, in ab- 
sence of magnetic order, cross the Fermi level, and are 
thus only partially occupied. For a system with partially 
filled ESs, in a ribbon, it might be energetically conve- 
nient to split those bands by inducing a magnetic order 
between the two ribbon-edges and thus opening an en- 
ergy gap, as in the well-known case of the zi edge. Fi- 
nally, we remind that the chemical reactivity of a given 
system can be quantitatively expressed as proportional 
to the electronic density of states near the Fermi level. 
As a consequence, all the structures having edge-states 
crossing the Fermi energy or in its immediate vicinity 
(independently on whether they are magnetic or not) are 
expected to be extremely reactive, since the edge-states 
are available to form chemical bonds. 

With these concepts in mind, we remark that the 
most stable structures at standard and high hydrogen- 
concentrations, namely an, Z211, and a22, have no edge 
states, in agreement with the experimental observation 
of large semiconducting gaps in narrow ribbons [6] . This 
implies that, besides being non-magnetic, they have a 
very low chemical reactivity. As a consequence, once 
formed, those edges are likely to be metastable well be- 
yond their respective windows of thermodynamical sta- 
bility in Fig. O On the other hand, magnetic edges, 
as the well studied zi, are stable only at extremely low 
hydrogen-concentrations and are also very reactive and 
thus unlikely to be stable at ambient conditions. 

We now show that the electronic structure of the edges 
can be interpreted in terms of aromaticity and accord- 
ing to Clar's rule, well known in organic chemistry since 
the 60's 2(3, 21 1 . The fundamental stability-criterion for 
hydrocarbons is that each carbon atom must have four 
saturated bonds. In organic chemistry, the so-called ben- 
zenoid aromatic ring can be defined as a resonance (du e 
to the derealization of tt electrons over the ring 2^, 23[) 
between 2 hexagonal rings with alternating single and 
double bonds as in Fig. [2]-E. For such a structure, all 
the bonds sticking out of the hexagon are single bonds 
and, as a consequence, two benzenoid rings cannot be 
adjacent. The energetic stability of a structure increases 
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with the number of possible resonances. As a conse- 
quence, since a structure containing n benzenoid rings 
displays 2" possible resonances, the most stable isomer 
of a given hydrocarbon is the one that maximizes the 
number of benzenoid rings (Clar's rule 0, [2l[). Bulk 
graphene is the ideal case of aromaticity since there are 
three equivalent representations in which all the tt elec- 
trons belong to a benzenoid ring and no double-bonds 
are present. In such a representation 1/3 of the graphene 
carbon-hexagons are benzenoid rings. 

The edge structures of Fig. [1] are displayed in a repre- 
sentation that maximizes the number of benzenoid rings. 
However, the resulting ring-density for some of the struc- 
tures (e.g. the Z21 or Z2111) is lower than the ideal 
graphene 1/3 density. We analyzed all the edges of Tab.H] 
and we indicate as "aromatic" those structures that have 
a density of benzenoid rings of 1/3 (e.g. the an, Z211 and 
8-22) ■ We note that all the aromatic structures do not have 
partially-filled edge-states (are non-metallic) and vice- 
versa. This empirical rule is valid for all the structures 
presently studied and can be understood in terms of the 
Clar's rule. In fact, whenever the ring-density near the 
edge is smaller than 1/3, there is a competition between 
the bulk, where aromaticity prevails with an exponential 
(2") weight, and the edge. It is precisely this competition 
that, imposing the ring density of the bulk, forces some 
edge carbon-atoms to have less or more than 4 saturated 
bonds and, thus, originates the electronic "defects", i.e. 
the edge states. The z 2 n and the z(600)nn ribbons 
are the only zigzag edges compatible with the aromatic- 
ity of the bulk. These edge configurations have in fact 
the same 1/3 periodicity along the zigzag edge as the 
benzenoid rings of "bulk" graphene. Moreover, mono- 
and di-hydrogenated armchair edges are both compatible 
with the 1/3 aromaticity of graphene, whereas asymmet- 
rically passivated edges, such as &21, are not. Notice that, 
for a given H density, the most stable edges are aromatic, 
if an aromatic structure is compatible with such a con- 
centration. The stability of fully-aromatic edges is con- 
firmed by the fact that polycyclic aromatic hydrocarbons 
whose structure can be fully represented by benzenoid 
rings (without double bonds) show particularly high sta- 
bility, high melting points, low reactivity [22J . and have, 
whenever possible, an or z(600)nn edges [2J, [24| • 

Concluding, the knowledge of the structure and sta- 
bility of the possible edges is a crucial issue to control 
the experimental conditions of the formation of graphene 
nanoribbons of desired properties. Here, we deter- 
mined the structure and stability of hydrogen-terminated 
nanoribbons. For reasonable hydrogen concentrations 
the most stable structures are the an, a 2 2 and z 2 n, which 
are not magnetic and do not present valence electronic 
states localized at the edge. These structures are, thus, 
expected to be non-reactive and metastable well beyond 
the region of their thermodynamical stability. In par- 
ticular, at high hydrogen concentrations graphene spon- 



taneously breaks into a 22 nanoribbons. Our results are 
rationalized by means of simple concepts of organic chem- 
istry that can be used to guide the search for the most 
stable edges in different chemical environments. 

We thank Th. Greber for discussions. Calculations 
were done at IDRIS (project n° 081202). 
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